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Abstract 

s * ■ The problem of a differential operator left- and right division is solved in terms 

of generalized Bell polinomials for nonabelian differential unitary ring . The defi- 
nition of the polinomials is made by means of recurrent relations. The expresions 
of classic Bell polinomils via generalized one is given. The conditions of an exact 
factorization possibility leads to the intertwine relation and results in some lineariz- 
ablegeneralized Burgers equation. An alternative proof of the Matveev theorem is 
CN | given and Darboux - Matveev transformations formula for coefficients follows from 

the intertwine relations and also expressed in the generalized Bell polinomials. 
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Problems of solitonics and, more general, of nonlinear partial differential equation 
(NLPDE) integration is simplifying when one introduce a special factorization of 



a (differential) linear operator associated with the problem. The technique was 
actively used in theoretical physics [1] as well as in mathematics from early [2] 
to [3] . First comes from attempts to build explicit solutions of the one-dimensional 
Schrddinger equation and leads to so-called " supersymmetry" and the second relates 
to general problems of a differential field theory. The "physical line" was developed 
for partial differential equations (PDE) in [4] where the direct connection with the 
theory of Darboux transformations was noticed. Both approaches connect factors 
with some Riccati equation that is a feature of the general factorization problem 
[5] of importance from algorithmic point of view [6]. The problem of factorization 
was succesfully investigated algebraically for a commutative field in [7,8] and non- 
abelian one in [9] that allows to find appropriate algorithm for computations [10]. 
Applications to nonlinear ordinary differential equations (NODE) decomposition 
theory shows a way to NLPDE via classical reductions of S. Lie (point symmetries). 
A search of approaches to solution of integrable NLDPE or looking for somewhat 
sort of a language in which algebraic algorithms could be formulated in a compact 
form lead us to the theory of Bell polinomials. The polinomials appears naturally 
in Darboux techniques [11] and even a disarmed eye see a special polinomial com- 
binations of derivatives in famous integrable equations. This way having in mind 
factorization of Schrddinger operators and Cole-Hopf linearization of the Burgers 
equation the Bell polinomials in fact appear in many contexts and generalized for 
nonabelian coefficients at [12]. The last result leads to a nonabelian Burgers hier- 
archy, directly linearized by the Cole-Hopf substitution analogue. From the side of 
integrable nonlinear differential equations form and their bilinear (Hirota) counter- 
parts the introduction of "binary" Bell polinomials allows to express the algorithm 
of linearization of a given NLDPE in compact combinatoric language [13]. 
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In the Section 2 of this paper we, after some notations, introduce left and right 
nonabelian polinomials by recurrence relations. The polimomials correspond to [12] 
but the form of their introduction is adjusted to our purposes and give some new 
useful relations. We also prove similar relationships in our context. In the Secton 3 
we make next step on a way of the generalization of bell polinomials. We introduce 
some auxiliary DO and show how to make a linear differential operator (LDO) 
division by left and right factors of the first order LDO and formulate conditions in 
which a factorization is possible (give a zero remainder) for a nonabelian field. Next 
section (4) links a solution of the problem of a factorization with generalized (and 
linearized by Cole-Hopf substitution analog) Riccati equation. At the last section 
we give a compact formulation and proof of the Matveev theorem [15,16] in our 
notations deriving the potenyials' (LDO coefficients) transforms in term sof Bell 
polinomials. 

2 Main notations and auxiliary results. Bell 
polinomials. 

Let K be a differential ring of the zero characteristics with unit e (i.e. unitary ring) 
and an involution denoted by a superscript *. The differentiation is denoted as D. 
The differentiation and the involution are in accordance with operations in K, i.e. 

1. (a*)* = a, (a + b)* = a* + b*; (aft)* = b*a*, a,b G K 

2. D{a + b) = Da + Db; D{ab) = (Da)b + aDb; 

3.(Da)* = -Da*; (1) 

4. the D n operators form a basis in K-module Diff(K) of differential operators. 
The subring of constants is Kq and a multiplicative group of elements is G. 

5. For any s G K there exists an element ip G K such that Dip = sip, that means 
also the existence of a solution of the equation 

D(/> = -<j)s (2) 

due to the involution properties. 

There are lot of applications in the theory of integrable nonlinear equations 
and quantum problems connected with rings of square matrices. In this case the 
matrices are parametrized by a variable x and D may be a derivative with respect 
to this variable or a combination of partial derivatives that satisfy the conditions 
(3,4). If D is such usual differentiation, then the involution "*" is the Hermitian 
conjugation. In the case of a commutator Da = [d,a], (Da)* = —[d*,a\. Having in 
mind such or similar applications we shall name the involution as conjugation. We 
do not restrict oneself by the matrix- valued case: an appropriate operator ring also 
good for our theory. 

Below we introduce left and right nonabelian Bell poplinomials (see also [12] 
The left differential Bell polinomials are defined by Definiton 1 

B (a) = e (3) 

and the recurrent relation 

B n {s) = £>B n _i(s) + B n _ 1 {s)s,n = 1,2, .. (4) 
Statement 1. If the element p G G satisfy the equation Dip = sp then 

D n p = B n {s)p, n = 0,l,2,.... (5) 
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Proof. Let us use the induction. At n = 0, 1 the equality (5) is trivial, from 
(3,4) n = 1 we get 

B 1 (s) = s; 

therefore, 

D ip = s <p = -Bi(s) ip. (6) 
Here and further we shall denote Dx = x', x € K for brevity especially when the 
operator D acts only on the nearest element. Let the equality (5) be trivial for some 
n, then, by means of the relations (6), (8) we obtain 



D n+1 ip = D(B n {s) ip) = B' n {s) ^ + B n {s) Dp> = B' n {s) ip + B n {s) s ip = 

(D B n (s) + B n (s) s) ip = B n+ i (s)ip; 

it means that (5) take place for the next value of n + 1, hence for all the natural 
n or zero. □ 

Evaluations by (6) give 

B 1 (s) = s; B 2 (s) = s 2 + Ds; B 3 (s) = s 3 + 2s' s + s Ds + D 2 s; 
B 4 (s) = s 4 + 3s' s 2 + 2s s' s + s 2 Ds + 3s" s + 3(Ds) 2 + s D 2 s + D 3 s. 



Right Bell polinomials are introduced similarily. 

Definition 2 The right nonabelian differential Bell polinomials are defined by 
the equalty 

B+(s) = e (7) 

and a recurrence 

B+(s) = -DB+_ l (s) + sB+_ l (s), n = l,2,.... (8) 

Statement 2. If the element <p G G satisfies the equation (2), then 

D n <j ) =(-l) n cj ) B+(s), n = 0,1,2,.... (9) 



Proof. Treat the induction again. At n = 0, 1 the equality (8) trivializes and 
n = 1 yields (2), whence (7), (8) give 

B+{s) = s; 

Otherwise, 

D4> = -(j)s = (f>B^(s). 



3 



Let the equality (9) is valid for some n, then, by means of (2), (8), (10) we get 

D n+1 <f> = D{(-l) n 4>B+{s)) = (-If (</>' B+(s) + <t)DB+{s)) = 
(-1)" (-0 s B+(s) + DB+(s)) = (-1)" +1 cP (-DB+(a) + s B+(s)) = 

so, it is valid for the next n + 1, hence for any natural or zero n. □ 
Calculations by (8) give 

Bf(s) = s; B£(s) = s 2 - Ds; B£(s) = s 3 - s' s - 2s Ds + D 2 s; 

= s 4 - s ' s 2 - 2s s' s - 3s 2 Ds + s" s + 3(Ds) 2 + 3s L> 2 s - D 3 s. 

The right polinomials go to ones determined at [12] if change D — ► — D . There 

exists a simple but useful link between the right and left Bell polinomials generated 
by the conjugation. 

Statement 3. The left and right Bell polinomials are connected by the following 
relations: 

B n (s)* = B+(s*), B+(s)* = B n (s*). (11) 

Proof. It is enough to prove the first one as the second is conjugate. For this 
proof the induction again starts from a trivial point n = of the equality (11 ), 
use the recurrences (4), (8) and give 

B n+1 ( s y = -DB n { s y + s * B n ( s y = -db+{ s *) + s * b+( s *) = s+ +1 ( s *). 

That means the validity for any n. □ 

If the ring is abelian, left and right polinomials coincide.. 

Remark 1. The statement 4 means, that for the Bell polinomials takes 
place a duality: any relation connecting right BP goes to one connecting the left 
BP and vice versa. Let us denote further 

L s = D-s. (12) 

Note that the recursion (7) may be written by means of the designation (12) as 

B+ +1 (s) = -L s B+(s), n = 0,l,2,...; 
with the simple corollary 

= (-!)" L?e, n = 0,1,2,.... (13) 
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3 One more generalization of Bell polinimials 



In the next section the problem of division of an arbitrary operator L by the 
operator L s will be studied. For the solution of it, we introduce auxiliary operators 
for a right division H n by means of 

Definition 3 

H = e. (14) 

and the recurrence:: 

H n = D H n -i + B n (s), n=l, 2, .... (15) 

Statement 5. The identity holds 

D n = H n ^L s + B n {s), n = l,2,.... (16) 

Proof. The induction gives the following. At n = 1 the equality (16) is trivial. 
Let it be for n, then 

D n+1 = D{H n _ 1 L s +B n {s)) = D H n -i L s +D B n (s) = D i? n _i L s +B n (s) D+DB n (s) = 

D L s + B n {s) L s + DB n (s) + B n (s) s = (D H n ^ + B n {s)) L s + B n+1 (s) = 

H n L s + B n+ i(s), 

it is (16), but after the change n for n + 1. Hence, the equality (16) is valid for all 
natural n. □ 

Coefficients of the operators H n are expressed via generalized Bell polinomials 
, that are defined by 

bf Definition 4 

B n ,o(s) = e, n = 0,1,2,..., (17) 

and recurrence relations 



B n ,k{s) = Bn-i,k(s) + £ ) #n-i,fc-i(s), k = l, n-1, n = 2,3, .... (18) 



B n , n {s) = DS n _i )n _i(s) + S n (s), n = l,2,.... (19) 
Statement 6. Generalized Bell polinomials are coefficients of the operators H n , 

i.e. 

n 

H n = Y, B n , n - k (s)D k , n = 0,1,2,.... (20) 

fc=0 
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Proof. When n = the equality (22) trivializes H = e. so that it is enough 
to establish for operators (20) the recurrence (15), by the following transformations: 



n n—1 

H n - D H n -i - B n (s) = ]T B n>n _ k (s) D k -J2D B n _i, n _i_ fc (a) D k - B n (s) = 

k=0 k=0 

n n—1 

E B n , n _ k {s)D k - E {B n _ 1>n ^ k {s)D k+1 + DB n ^ n ^ k {s)D k )-B n {s) = 

k=0 k=0 

n n—1 n—1 

E B n , n - k ( 8 )D k - E B n -i,n-i-k(s) D k+1 - E DB n _ hn ^ k (s)D k - B n (s) = 

k=0 k=0 k=0 

n n n—1 

E B n , n _ k (s)D k -J2 B n _ hn _ k (s)D k -Y DB n - 1 , n _ 1 _ k (s)D k -B n (s) = 

k=0 k=l k=0 

n-1 

(B n , (a) -B n _ h0 )D n + E (Bn,„-fe(s) - B B -i, re - fe (s) - M n -l,n-l-fc(s)) ^+ 

k=l 

{B n , n (s) - £>B n _i, n _i(s) - B n (s)) = (e-e)D n = 0. 



As the equality (14) and the recurrence (15) define the operators H n , in unique 
way, the result of the calculations means the validity of (20). □ 

The formulas (18) and (19) are simple, but nor very convenient for evaluations 

of B njk (s) 

(you forced to go "ladder" way) , therefore we suggest more complicated but prac- 
tically easier algorithm. For this we put the representation (20) into (18); getting 



= D n+1 -H n L s -B n+1 (s) = D n+1 -J2 B n ,n-k(s)D k (D-s)-B n+1 (s) = 

k=0 

n n 
D n+1 -Y / Bn,n^(s)D k+1 + Y / Bn,n-k( S )D k S-B n+1 ( S ) = 
k=0 k=0 
n+1 n k /,\ 

D n+i _ £ B n>n _ fc+1 ( a )£>* + E B n ^ k {s) £ . D^sD* - B n+1 (s) = 

k=l k=0 i=0 V / 



(e-B n , (s))D n+1 -Y B n ^ k+1 {s)D k + ]T ( ) ) B n , n _ k {s) D k - l s D l -B n+1 {s) 

k=l 0<i<k<n 



I 



(e-B n , (s))D n+1 -Y B n , n _ k+1 (s)D k + E ) B n , D k -B n+1 (s) 

k=l < k < i < n * ' 



n n n / ■ 

I 



( e _ e ) _ ^ B n , n _ fc+1 ( S ) £>* + E £ J Sn.n-iCs) Z^Z^ - B n+ i(*) 

k=l k=0 i=k V ' 



E -5n,n-fc+l(s)+E [I.)5„ 1 „-i(a)D < - fc s] D fc +(E ^n,n-i(s)^S-B„ +1 (i 
k=l \ i=k W ) \i=0 



The following formulas are extracted: 



£n,n-fc+i(s) = J2 { I) B^ n ^{s)D l ~ k s, k = l,n, n = 0,1,2, ... ; (21) 



i=k 



B n+1 {s)=Y J B n ^ l {s)D i s, n = 0,1,2,.... (22) 

i=0 

The formulae (22) expresses the standard (nonabelian) Bell polinomials via the 
generalized ones; it may be rewrited as 



B n+1 (s) =J2 B n ^{s)D n - l s, n = 0,1,2,.... (22) 



If in the formulae (22) rearrange the summation by k — > n — k + 1, then after 
simple work yields: 



fc-i 



n — % 



B n,k(s) = 2^ y n _ k + ij B n,i(s)D K ~ l ~ l s, k = l,n, n = 0,l,2,.... (23) 
Evaluation of the generalized Bell polinomials by (22) gives 



B n:1 {s)=s; B ni2 (s) = s 2 + nDs; 5 n , 3 (s) = s 3 + n s a + (n-l) s Ds+ I I D 2 s; 



B nA (s) = s A + n s' s 2 + {n - 1) s s' s + (n - 2) s 2 Ds + ^ J s" s + n(n - 2) (Ds) 2 + 



For a solution of problem of the left division of a differential operator, L by L s 
the similar but simpler consideration is necessary. The analog of the Statement 
5 is 

Statement 7. The following identity is valid 

D n = L s H+_ l + B+{s), n = l,2,..., (27) 

where 



H+ = ^2 B n-k(s)D k , n = 0,1,2,.... (28) 

k=0 

Proof. Let us solve the problem of the left division of the operator D n by the 
operator L s i.e. put 



D n = L s H+_ 1 + r, 

where r -remainder, and 



(29) 



n-l 

H+_ 1 = J2hD k . (30) 

fc=0 

The remainder, r and the coefficients b k are to be determined. 
Substituting (30) into (29) and, transforming, we arrive at 



n— 1 n— 1 n— 1 

L s i/-_ 1 +r = ( J D- S ) ^ 5 fc D fc + r = £ Db k D k sb k D k + r 

k=0 k=0 k=0 

n—1 n—1 n—1 n—1 n—l 



^2(b k D k+1 +Db k D k )-Y, sb k D k +r = Y J b k D k+1 +Y J Db k D k -^ sb k D k +r = 

k=0 k=0 k=0 k=0 k=0 

n n—1 n—1 

]T D k +Y / (Db k -s b k ) D k +r = (Db -s 6 +r)+^ (& fc _i+£%-s b k ) D k +b n _ 1 D n 



k=l k=0 k=l 

n-l 



(L s b + r) + J2 (h-i + L s b k ) D k + 6 n _x D n . 



k=l 



Comparing the result with the left-hand side of (31), we obtain 



6 n _i = e; b k -i = -L s b k , k = 1, n - 1; r = -L s b . 
Recalling (7) and the recurrence (8), we represents: 

&fc = -B;U-i( a )> fc = 07¥^T; r = B+(s). 
The equalities (24) and (25) are their obvious corollaries. □ 

3. Division and factorization of differential operators. 
Generalized Riccati equations. 

Let 

N 

L = Y / a n D n , a n eK,- (28) 

?1=0 

be a differential operator of the order N. We shall study a right and left division 
of L by the operator L s „ defined by the equality (12), Suppose 

L = ML s + r, (29) 

L = L s M + + r + , (30) 
where M, M + -results of right and left division and r, r + - remainders. 
The statement 5 allows to solve the problem in a simple way. 
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Statement 8. If the representation (29), is valid, then for the remainder r 
and the result M yields: 

N 

r = Y1 a n B n (s); (31) 

ra=0 
N 

M=^a n H n . 1 , (32) 

n=\ 

or 

N-l 

M=Y,bnD n , (33) 

n=0 

where 

N 

b n = E fl * B *-i,nW. n = 0, JV-1. (34) 

fc=n+l 

Proof. Multiplying both sides of the equality (16) by a n and summing by n, one 
get (29), (31) - (33). After the substitution of the representation (22) into (34) it 
results 

N N n-l 

M= = °™ S B n-l,k(s)D k = ^ O n k ( S ) D k = 

n=l n=l k=0 l<n<N,0<k<n-l 

N-l N 

E a k B k _ hn (s)D n =J2 E afc£ fc _ lin ( S )£> n , 

l<k<N : 0<n<k-l n=0 fc=n+l 

from what the formulas (33), (34) follow. 
As a corollary one get 

Statement 9. For the linear operator L to be a right divisible by L s without 
remainder it is necessary and enough, that s to be a solution of the differential 
equation 

N 

Y / a n B n (s) = 0. (35) 

n=0 

// this condition holds, the operator L factors as : 

L = ML S , 

where the value of the operator M is given by (32) or by the expressions (33), 
(34). 

The equation (35) is nonlinear. At N = 2 it is Riccati equation, therefore it 
is natural to name it as a generalized right Riccati equation. The right Riccati 
equation is generalized by means of the Statement 2. 

Statement 10. Let an invertible function if be a solution to the linear differ- 
ential equation 
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N 



E an D n v? = 0. (36) 



n=0 



Then the operator L, defined by the equality (31), is right divisible by L s , where 



s = If' ip 1 . 



Morover s is a solution of the right Riccati equation (35). 

For the solution of a left-division problem letus write a result in a form 

JV-l 

M + =Y,KD n , (37) 

?1=0 



looking for the determination of &+, n = 0, n — 1. 

Substitute the representation (37) into the right-hand side of the equation (30). 
By means the calculus analogues to that was used for the proof of the Statement 7, 
you obtain: 

= a N ; (38) 
b+ = a n+1 - L s b+ +1 , n = 0,N-2; (39) 
r+ = a -L s b . (40) 

Solving subsequently the equations (37) - (40) one goes to 

N 

bt= E (-l) fc - n_1 ^- n - 1 a fc , n = 0, JV-l; (41) 



fc=n+l 



AT 

r+ = ^(-l) fe L^a fe . (42) 

fc=0 



The entities b£, n = 0, N — 1; r + may be expressed in terms of thr right Bell 
polinomials if use the equality (12) and take into account 

L k s a = L k s ea={-l) k B+(s)a. 
Hence links (41), (42) transform to: 

N 

bn~= E £ fc + -„-i(*K> n = 0,iV-l; (43) 

k=n+l 

N 

r + = J2 B k(s)ak- (44) 

k=0 

Formulas (30), (37), (43), (44) (or (41), (42)) give a solution of the left division 
problem of L by L s . So , it is proved 
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Statement 11. if the representation (31), is valid, then for the reminder r + 
and the result M + the formulas (44) (or (42)) and (37), take place and for the 
coefficients of the operator M + there is a representation (43) (or (41)). 

The straight corollary of this sentence is the following 

Statement 12. For the operator L to be left divisible by the operator L s 
(without remainder) it is necessary and enough, that s be a solution of the differ- 
ential equation 

N 

B k (s) + a k = 0. (45) 
if this condition holds, then the operator L factorizes as: 

L = L s M + , 

where the value of the operator M + is given by the expressions (37), (43) (or 
(42)). 

The nonlinear equation (45) is called as generalized left Riccati equation. 
The left Riccati equation is linearized obviously by the Statement 3. As a result 
we have 

Statement 13. Let an invertible function (p satisfy to the linear differential 
equation 

N 

Y / (-l) n B n (s) + a n D n ^ = 0. (46) 

n=0 

Then the operator L, determined by the equality (31), is left divisible by the operator 
L s , where 

s = —tp^ 1 ip' . 

The function s is a solution to the left Riccati equation (45). 

Remark 2. Following F. Calogero classification, the Riccati equations (38), 
(45) are C-integrable dynamical systems. 

4. Darboux transformation. Generalized Burgers equations. 

We shall show that the problem under consideration of the operator division is 
connected directly with Darboux transformation. For the purpose we take a version 
in which at the ring K there exists one more differentiation Dq, which commute with 
the operator D D, 

D D = DD . 

It may be a differentiation by a parameter t. 
Let us list auxiliary commutation relation. 

L s r = rL s + Dr+[r,s]. (47) 
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Really, 



L s r—rL s = {D—s) r—r {D—s) = D r—sr—r D+rs = r D+Dr—sr—r D+rs = Dr+[r, s] 

Taking into account the equalities (48) and (31), we calculate 

L s (D -L) = L s Dq-L s L = D L s +D s-L s (M L s +r) = D L s +D s-L s M L s -L s r = 
D L S + D s-L S M L s -rL s -Dr-[r, s] = {D - L s M -r) L s + D s - Dr - [r, a]. 

and arrive to the following relation: 

L s {D — L) = {D -L)L S + D s-Dr- [r, s] , (48) 

where 

L = L s M + r. (49) 

We results in the important conclusion: 

Statement 14. // a function s satisfy the equation 

D Q s = Dr + [r, a], (50) 
the operator L s intertwine the operators Do — L 3 Dq — L, 

L S {D -L) = {D -L)L S . (54) 
Now we would obtain the explicit expression for L,in the terms of (49), (31), (32); 
namely 

N N 

L = L s M + r = {D-s)J2 a n H n -i + a ™ B ^ = 

n=l n=0 
N N N 



Da n H n ^i -J2sa n H n _i + ^ a n B n {s) = 

n=l n=l n=0 

N 

^ {Da n H n _i + a n D H n _i - sa n H n _i + a n B n {s)) + a B {s) = 

n=l 

N 

{Da n H n _ x + a n {H n - B n {s) + B n {s)) - s a n F n -i) + a e = 

n=l 

N 

^ ( a n Hn-i + a n H n — s a n i? n -i) + ao- 

n=l 

Finally, the transfomed operator is 

12 



N 



L = ^2(a' n R n -\ + a n H n - sa n i? n -i)) + ao- 



(52) 



n=l 



Let us write the equation (50) in the explicit form by virtue of the formula (31). 



N 



N 



a n B n (s), s 

.n=0 



Dr + [r, s]=DJ2 a n B n (s) + 

n=0 

N N 

^2 Da n B n (s) + Y [a n B n (s), s) = 

n=0 n=0 
N 

Y (a n B n(s) + a n DB n (s)) + ^ (°« B n(s) s- sa n B n (s)) 

) n=0 
N 

Y (a n B n (s) + a n DB n {s) + a n B n (s) s - s a n B n {s)) = 

n=0 
N 

E {a'n B n{s) + a n (DB n (s) + B n (s) s) - sa n B n (s)) = 



N 



n=0 



n=0 



N 



Y ( a 'n B n(s) + a n B n+1 (s) - sa n B n {s)). 

n=0 

It is established that for the intertwine relation (51) validity it is necessary and 
enough that s should be a solution of the equation 



N 

E 

n=0 



D s = V (Da n B n (s) + a n B n+ i(s) - sa n B n (s)). 



(53) 



Remark 3. The equation (53) is nonlinear, but linearizable, i.e. (by Calogero 
classification) it is C-integrable. This equation (in different forms) was introduced 
at [12,16]. The form we suggest is most compact nad convenient for a further 
investigations, e.g. in the franework of bilineraization technique of [13]. 

Example. 

Let 



L = D\ 



Then 



L = H 2 -sH l = B 2 ,o{s) D 2 + B 2 ,i(s) D + B 2 , 2 (s) - s (Bi,o(s) D + B lt 1 (a)) 
D 2 + sD + {s 2 + 2Ds) -s(D + s) = D 2 + 2Ds. 

The equation (53) takes the form: 



D s = D 2 s + 2Dss. 
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In the case of scalar functions it is known Burgers equation. By this reason and due 
to the eqaution are integrable by the Cole-Hopf transformation, the equation (53) 
is naturally named generalized Burgers equation. 

Statement 15. Let an invertible function ip, is a solution to linear differential 
equation 

D <p = Ltp. (54) 
Then the function s satisfy the generalized Burgers equation (56). 
Proof. Let us note that 

(Dr +[r, s])<p= (L s ) r ip + r sip = (L s ) r p> + r D tp = L s (r ip) = L s Lip>. 

Further, acting to the equation (53) by the operator L s and accounting for the 
relations (45) and (2), (16), we have 

= L s (D (p - Lip) = L s D ip - L s L(p = (D L s + D s) ip - L s L(p = 

D L s ip + D sip - (Dr+ [r, s])ip = (D s - Dr - [r, s])ip, 

due to the existence of the inverse element for <p one obtain (50), that is equivalent 
to (52). □ 

The obvious corollary of the intertwine relation (51) and the statementl6 is 
Theorem (Matveev). Let functions ip and V are solutions of the equations 

D Q ^ = Lip, D ip = Lip, 
for an invertible <p Then the function 

ip = L s \p = Dip — stjj, s = Dip>ip~ 1 , (55) 
is a solution of the equation 

D ip = Lip. (56) 

The last statement accomplishes the proof of Matveev theorem [15,16]. The 
equality (55) gives a representation of the transformed operator in terms of the 
generalized Bell polinomials. The explicit expression for the transformed coefficient 
is 

a N [l] = a N , 

N 

a k [l] =a k + Y^ [a n B n ,n-k + (<4 ~ S«n) 5 n-l,n-l-fe], k = 0, N - 1. 
n=k 

5. Conclusion . 

It is shown that the division procedure for linear differential operators naturaly 
leads to the solution of its factorisation problem thats links to intertwine relations 
and Darboux transformations. The representation constructed here may, perhaps, 
open new possibilities to build up and study new integrable systems. One of us 
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Appendix. Evaluation of generalized Bell polinomials.. 

Calculations by the relation (26) give the following: 
At n = 1 

BnM = E ("J"/) B nA s ) (D-'s) = ( n ) B n , (s) (D°s) = s; 



i=Q 



At n = 2 



Bn,2{s) = jz ( n i _ l )B n , i (s)(D 1 - i s)= (f\ B n , (s)(D i s)+( n Q ^ B n , 1 (s) (D°s) 



n Ds + s s = s 2 + n Ds. 



At n = 3 



2 / A 

n — % 



i=0 



Q B n , (s)D 2 s + 1 ^ B n , 1 (s)D 1 s+ 2 j B n>2 (s)D°s = 
U j D 2 s + (n - 1) s Ds + (s 2 + n Ds) s = s 3 + n Ds s + (n - 1) s Ds + I ™ \ D 2 s. 



At n = 4 



3 / A 
n — i 



B nAs) = E [l_ i )B n , i (s)(D s - i s) = 



™jB nt0 (s)D 3 s+(^ 2 1 jB ntl (s)D 2 s+(^ i 2 jB nt2 (s)D l s+ 

( n 3 )^n,3( S )^ = 



|^s 3 + n Ds s + (n - 1) s Ds + ^ J L> 2 sJ s = 
™ J D 3 s+ Y s £> 2 s + (n - 2) s 2 + n(n - 2) (Ds) 2 + S 4 + tifls s 2 + 
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( 



n — 1 



2 



) 



sD 2 s+ 
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